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Abstract. We study energy level statistics of a critical random matrix ensemble of 
power-law banded complex Hermitian matrices. We compute the level compressibility 
via the level-number variance and compare it with the analytical formula for the exactly 
solvable model of Moshe, Neuberger and Shapiro. 



1. Introduction 

Random matrix ensembles (RME) provide a natural frame work for the statistical 
description of quantized energy levels of complex systems. This is probably best 
explained by paraphrasing Dyson [T: in ordinary statistical mechanics a renunciation 
of knowledge about the system is made i.e., by assuming that all states of a large 
ensemble are equally probable. Therefore only the overall behavior of the system is 
accessible. Thus orthodox statistical mechanics is clearly inadequate for the discussion 
of energy level statistics of complex systems. In energy level statistics one wishes to 
make statements about the fine detail of the energy level structure, and such statements 
cannot be made in terms of an ensemble of states. Random matrix theory (RMT) turns 
out to be a suitable tool in making such statements about the eigen-spectrum of complex 
systems e.g., energy levels of highly excited nuclei [2|. The underlying assumption in 
RMT is that one renounces exact knowledge not of the state but of the nature of the 
system itself. This is the fundamental difference between ordinary statistical mechanics 
and energy level statistics. 

There are three RME that emerge from an analysis of real eigen-spectra of random 
Hermitian matrices 0; the Gaussian orthogonal ensemble (GOE) for a system with 
time-reversal without spin-| interactions, the Gaussian symplectic ensemble (GSE) 
for a system with time- reversal and spin-| interaction, and the Gaussian unitary 
ensemble (GUE) for a system without time-reversal. These ensembles (generally called 
Wigner-Dyson (WD) statistics) exhibit strong energy-level correlations leading to the 
phenomenon of level repulsion ^ in their respective eigen-spectra. One could add the 
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Poisson statistics (PS) as the fourth ensemble which does not exhibit level repulsion in 
it's corresponding eigen-spectra. 

In general, the WD statistics describe a system of delocalized states and the 
PS describe a system of localized states. The transition between these two cases 
controlled by the strength of disorder [i.e the Anderson model [SI) is a quantum critical 
phenomenon. At a critical disorder the states acquire a special property of multifractality 
inilZ] which marks a qualitative difference as compared to truly delocalized states in a 
metal and truly localized states in an insulator. These critical states correspond to the 
critical level statistics (CLS) [311] • 

There are numerous suggestions in literature P ITI H ITT ] IT ^ ITH | I7t IT U ITK | ITH | ITT] 
describing CLS by ensembles of random matrices. Initially, there were two main 
requirements for the critical RME: the first was that the level statistics be system- 
size independent and secondly, the ensemble should be able to interpolate between WD 
and PS |H]. Later on it was realized [UlEl that the key property of the critical states 
relevant for the CLS is multifractality. The power-law banded RME suggested in [TT] 
emerged as the best candidate for the critical RME as its level statistics is known to 
possess these two properties and its eigenfunction statistics is multifractal. Thus it has 
been suggested to be the critical random matrix ensemble (CrRME). In the limit 
of large bandwidth this model turns out to be the best deformation of the WD-RME 
that generates correlations in the eigenstates but retains all the basic properties of WD 
theory, including the factorizability of higher spectral correlation functions and mapping 
onto the Calogero- Sutherland model of fictitious interacting fermions 

In the absence of an exact solution of CrRME various mappings on and approximate 
correspondences with exactly solvable models have been used in order to study the 
properties and CrRME. Here we mention the mapping onto the nonlinear sigma- 
model PT], and the approximate equivalence jH] between the CrRME and the exactly 
solvable RME suggested in Refs. jTSHEj. Making use of the exact solution [12] it was 
possible j2ni to compute the level spacing distribution function (LSDF) for this model 
and compare it with the exact diagonalization of the 3d Anderson model for the critical 
value of disorder. The coincidence of results (for certain values of the parameter b which 
were found from the fitting of the far tails of LSDF) was amazing |2n]l23j : in all three 
Dyson symmetry classes the deviation was not larger than between the Wigner surmise 
and the exact WD distribution function. 

In this paper we study the level number variance in the CrRMT suggested 
The ultimate goal here is to compare the numerical simulations of this critical ensemble 
and the analytical results for the exactly solvable ensemble in Ref. 110] (which will be 
referred to as MNS) for all values of the control parameter b. For large 6^1 (weak 
multifractality) these two ensembles are known to be equivalent jH]. The first correction 
to the Poisson spectral statistics at small 6 <^ 1 (strong multifractality) is also the 
same pT] l22j. We would like to check these predictions numerically for the unitary 
CrRME and to see how large are the deviations between these two RME for 6 ~ 1. 

The paper is organized as follows: in Section |21 we give a general overview of the 
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level-number variance (E2-statistics) and the calculation of the level compressibility, in 
Section El we give some numerical details and finally a summary and outlook of this 
work is given in Section 0] 



2. S2-Statistics and Level Compressibility 

Our model is the power-law random banded matrix defined as an ensemble of random 
Hermitian N x N matrices belonging to GUE. The real and imaginary parts of the 
matrix elements are independently distributed Gaussian random variables with zero 
mean (Hij) = and the variance of the matrix elements is given by; 

where 6 is a parameter that controls the fractal dimensions of the critical eigenfunctions 
and /3 is the symmetry parameter i.e., /? = 1,2 and 4 for the GOE, GUE and GSE 
respectively. This model remains critical for any arbitrary value of b. 

From a physical point of view the GUE can be realized by putting a system in a 
strong external magnetic field i.e., provided the splitting of levels by the magnetic field 
is of the same order of magnitude as the average level spacing in the absence of the field. 

In this paper we study the level compressibility b) via the level- number 
variance S2 {{n)) of the above model. The level-number variance is a statistical quantity 
that provides a quantitative measure of the long-range "rigidity" of the energy spectrum 

0; 

E,{{n)) = {n')-{n)' (2) 

The behavior of S2 is known (as a function of the mean number of levels (n) ^ 1) in 
the delocalized and localized phases respectively. However, in the critical regime, the 
level-number variance has been conjectured to be Poisson-like PUCB]; 

^ In ( (n) ) delocalized, 
X {n) critical, (3) 
{n) localized, 

where the level compressibility x is another important parameter to characterize the 
localization-delocalization transition (LDT) and takes values < x ^ I5 being zero in 
the delocalized state and unity in the localized state. Formally, one computes x as; 

X ~ lim lim — ;\ " 4 

(„)^oo Af^oo d{n) 

The limits in Eqn.(0]) do not commute and their non-commutativity is attributed to 
the violation of the normalization sum rule ^H]- It is remarkable (see [23111011122] for 

details) that x can be expressed in terms of the fractal or the correlation dimension 
D2 and the spatial dimension d. In Ref. |18| an explicit and very simple relationship 
between the critical spectral and eigenfunction statistics has been suggested: 

x4(i-^) («) 
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There is vast analytical and numerical evidence that this formula is correct at weak 
multifractality. However, its status for the case of strong multifractality (small D2) 
raises much doubts We later see that whilst our numerical data are in agreement 
with Eqn.(j21) in the limit they are in disagreement pT|l^ for the opposite limit. 

At the same time the overall agreement with the analytical result Eqn.(lO) for the MNS 
model is unexpectedly good for all values of h and the moderately large matrix sizes up 
to = 400. We were not able to detect any systematic dependence on the size of the 
matrix in this range. 

In the next section we turn to the numerical details of how to compute the level 
compressibility. 



3. Results and Discussion 



We performed numerical simulations of the RME described by Eqn.(^ for values of 
h^^ G [0.4,40]. The matrices are not sparse and as such the use of fast and efficient 
algorithms i.e., the Lanczos algorithm is not possible. However, we calculated the 
spectra of N x N matrices for A^ = 50 with 10000 realizations to A^ = 400 with 
2000 realizations by exact numerical diagonalization and hence accumulating ~ 1 x 10^ 
eigenvalues for each value of A^. The results we discuss here are for the GUE only. 
The average density of states as a function of energy E is well described by; 

'~^bj, l + (l/z)expP/2) (E^ + t^)] 
where z = exp {it (3b) — 1. The two limiting cases of Eqn.(jni) are; 

4i^2T(b - E^ Semi-circle if irBb > 1 

piE) = < : rj_ , ^ . (7) 



A^y^exp {—l3E'^/2) Gaussian otherwise 



All our calculations are done at the band center (containing (n) ^ 1) where the density 
of states is nearly constant. The calculation of x from S2 is non-trivial as there are 
finite-size effects ~ — (n)^ /N in S2 originating from the long tail R{s) ~ ~1/N of the 
two-level correlation function. This is circumvented by fitting S2 data with a polynomial 
in (n) of degree m: 

m 
k=2 

Then we identify the coefficient of the linear teTm as x and xo is the constant teim. Both 
these terms are expected to have only weak A^-dependence. This is in contrast to that of 
the coefficients Xk oc 1/N'^~^ with A; > 1 that have a strong A^-dependence which describe 
finite-size effects in S2. By making the polynomial fit Eqn.(jH)) we separate the finite-size 
effects in E2 and the weak A^-dependence in x and in xo- This weak A^-dependence in 
X and in xo could be of fundamental interest and needs further investigation. 

In Fig^we plot S2 for A^ = 400 and S2 {{n)) is linear as a function of (n) as is 
expected for a critical ensemble. The lines through the data points are in agreement 
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Figure 1. The typical behavior of the level- number variance S2 {{n)) for N — 400 
and for b^^ = 0.4, 1.2, 2.0, 2.8 and 3.6. The lines through the data points correspond 
to Xo {b) + (j^) {n). The value of xo is on average N independent (see Fig|3Il. The 
linear behavior of S2 is in perfect agreement with the theoretical prediction of (2111 18| 
(see Eqn.®). For clarity in the numerical data we skip every Srrf symbol. 



with the Poisson-hke behavior i.e., S2 ~ x {n) + xo- For large h ^ l/27r the value of 
X is in good agreement with Eqn.(jn)) and Eqn.Q that follows from Eqn.dS)) and the 
correlation dimension D2 found in Refs. f3J|2Hl '■ 

X = ll^T^h (9) 

However, our results are in disagreement with Eqn.(0) for small h <^ l/27r pTll22j . The 
analytic formula for x (MNS in FigI2I) that follows from JU] for all values of h reads: 

^ ^ln(J) ^ PolyLog [-1/2, 1 - exp (27rb)] 
^ d\Yi{z) PolyLog [+1/2, 1 - exp (27rb)] ^ ' 

where 

dp 

(1 + exp 

and PolyLog [n,z] is an analytic continuation of YlT=i^^/^"^- This function can easily 
be computed by use of mathematica. 

The main result of this work is that the values of x (for modestly large < 400) 
of the critical ensemble due to |TT] agree unexpectedly well with the theoretical formula 
due to jmi, Eqn. (fTn|l . This agreement is, however, not perfect since the data points lie 
below Eqn. ()10|) in the range 3 < b^^ < 10, where the function x (^~^) Eqn.()10|) has a 
peculiar curved shape. 





Figure 2. The spectral level compressibility x as a function of b^^ and N. The 
values of x were calculated by a fit-polynomial of degree m = 3. The solid line is the 
analytic formula (MNS) due to jTH] and the dashed line is Equ.®. x is on average iV 
independent. 




Figure 3. The constant term xa{N,b) of the fit-polynomial with degree m = 3 in 
Eqn.®. Xo decays to zero for values of < 10. The data here is the same as in 
FigEl 
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4. Conclusion 

We have been able to compute x as a function of b^^ for the GUE critical ensemble. 
Our results are in agreement with the theoretical formula Eqn. pOj) except for the range 
3 < b^^ < 10 where the numerical data are below Eqn. ()10|l . 

Furthermore, xo and x are almost A^-independent which is a signature of criticality. 
It is important to note that our simulations exclude the strong A^-dependence of x but 
they do not exclude the weak ln(A^)-dependence of x which is not trivial to realize 
numerically. Recently, it has been shown analytically |2Z] that in the CrRME the 
corrections to x of order 6^ have a logarithmic A^-term i.e —b'^\n{N) in the crossover 
between the GOE and the GUE. In both the GUE and the GOE this logarithmic 
factor is absent up to order b^. It's presence is not excluded in higher b corrections, 
e.g as — 6^1n(A^). Such a logarithmic dependence is of fundamental interest. This is 
because the presence of ln(A^)-corrections would signal the incomplete criticality of the 
model. The logarithmic terms may drive the system to a different, stable, set of critical 
points where the x niay be drastically different from the MNS result. Unfortunately, the 
accuracy of our simulations is not enough to detect such a weak logarithmic dependence. 
This point needs further investigation. 

We note also that the constant term xo carries important information about CLS 
since it has a finite limit as — > cx). As is seen in figure El the value of xo as a function 
of b~^ decays to zero with increasing b^^ which is (so far unnoticed) an important 
characteristic feature of CLS. 

In summary we have been able to compare via numerical simulations of the GUE 
critical ensemble ^T] and the theoretical formula Eqn. pO|) . The two results agree but 
in the region 3 < b^^ < 10, the simulations are below the theoretical calculations. 
However, we emphasize that our results are preliminary and further work is in progress. 
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